subroutine comput_flux1
use var
implicit none

    integer m,j
    
        do m=1,4
            call flux1_m(m)
            call flux1_p(m)
            do j=0,gridsx
                flux1(j,m)=fp1(j,m)+fm1(j,m)
            end do            
        end do


end subroutine